Low energy electronic states and triplet pairing in layered cobaltates 



in 
o 
o 

(N 

o 
O 

m 



X3 
O 

o 



in 
o 
o 

in 
o 



o 



X 



G. Khaliullin,!'^ W. Koshibae,^ and S. Maekawa^ 

^ Max-Planck-Institut f-iir Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 
E. K. Zavoisky Physical-Technical Institute of the Russian Academy of Sciences, 420029 Kazan, Russia 
Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 
(Dated: February 2, 2008) 

The structure of the low-energy electronic states in layered cobaltates is considered starting from 
the Mott insulating limit. We argue that the coherent part of the wave-functions and the Fermi- 
surface topology at low doping are strongly influenced by spin-orbit coupling of the correlated 
electrons on the t2g level. An effective t-J model based on mixed spin-orbital states is radically 
different from that for the cuprates, and supports unconventional, pseudospin-triplet pairing. 
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The layered cobalt oxides Naa;Co02 exhibit a number 
of remarkable properties with potential applications. At 
large x, they are unconventional metals showing a large 
thermopower (suppressed by magnetic field 0) and also 
signatures of localized magnetic states. Such a mixture of 
itinerant transport and Curie- Weiss magnetism ("Curie- 
Weiss metal" 0]) is not easy to reconcile, and there 
is growing evidence that special charge/orbital ordering 
correlations are at work at large x Static charge or- 
dering is observed at x = 0.5 01 • These observations 
are not totally unexpected, as a rich interplay between 
spin/ charge /orbital degrees of freedom is a common phe- 
nomenon in transition metal oxides. Yet a real surprise 
and new challenge for theory was the recent discovery jj| 
of superconductivity at Tc ^ 5 K in water-intercalated 
cobaltates with low sodium content about x ~ 0.3. It 
is believed that the superconductivity of the cobaltates 
emerges from a correlated metallic state with enhanced 
electronic mass and it may have an unconventional, 

possibly spin-triplet 0, iS, pairing symmetry. 

Theoretically, regarding Na2;Co02 as spin 1/2 Mott 
insulator doped by spinless charge carriers, a t-J model 
similar to that for the high- Tc cuprates has been consid- 
ered However, the predicted time rever- 
sal violatingdi_-|- id2 singlet state is not supported by 
experiment [6|, |9j , and other proposals employing fermi- 
surface nesting and /or charge-density fluctuations have 
been discussed 0, 0| in favor of spin-triplet pairing. 

In this Letter, we show that the relevant t-J model 
for the cobaltates is in fact qualitatively different from 
its simplest version used in Refs. 
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key point is that the low-energy electronic states of the 
C0O2 layer are derived from the Kramers pseudospin 
doublets of the Co^^ ion with mixed spin and orbital 
quantum numbers. The projection of the (initially spin- 
antiferromagnetic) J-interaction on the pseudospin states 
has nontrivial consequences for the internal structure 
of the Cooper-pairs, leading to novel pseudospin-triplet 
pairing. 

The basic structural element of Naa;Co02 are C0O2 
layers, which consist of edge sharing CoOg octahedra 



slightly compressed along the trigonal axis The Co 
ions form a 2D triangular lattice, sandwiched by oxygen 
layers. Sodium doping is believed to introduce spinless 
Co'^"'" states into the spin 1/2 Co^+ background. 

A minimal model for the cobaltates should include the 
orbital degeneracy of the Co'^^-ion E3i where a hole in 
the d^{t2g)-she\\ has the freedom to occupy one out of 
three orbitals a = dyz, b = dxz, c = d^y- The degeneracy 
is partially lifted by trigonal distortion, which stabilizes 
the Aig electronic state {a+b-\-c) / \/3 over the iJ^-doublet 
{e^^'^a -f e=F«"^& c)/V3 (hereafter ^ = 27r/3): Ha = 
A[n(E') - 2n{Aig)]/3. The value of A is not known; 
Ref.[l3| estimates it ~25 meV, while the band structure 
calculations for a related structure give ~100 meV E3- 

In terms of the effective angular momentum leff = l of 
the functions Aig and E'^ correspond to the 
\lz = 0) and \lz = ±1) states, respectively. Therefore, a 
hole residing on the E'^ orbital doublet will experience 

an unquenched spin-orbit interaction Hx = —X{1 ■ s). 
The constant A for Co'*+ is about 640 cm^^ « 80meV 
|20| . Although A is somewhat smaller than the bare hop- 
ping matrix element t OA eV in cobaltates (inferred 
from the free-electron bandwidth ~1 eV the spin- 

orbit coupling strongly affects the coherent motion of the 
fermions when the quasiparticle bandwidth is reduced by 
correlation effects to much smaller values ~ 0.1 eV li^. 

At X — 0, H — Ha + H\ is diagonalized by a transfor- 
mation: 
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ICg 



(1) 



where cg — cos 9, sg — sin 6, a — (a, 6, c) and tpa = 
(0, — (/), 0), correspondingly. The angle 9 is determined 
from tan26' — 2-\/2A/(A -I- 2A). As a result, one obtains 
three levels, fg,ga,hg [see Fig^^a)], each of them are 
Kramers doublets with pseudospin one-half a. The high- 
est, /-level, which accommodates a hole in the t2g-shell, 
is separated from the f;-level by Sf — Sg — A-t- i(A/2-|- 
A)(l/cos26' - 1). This splitting is ~ 3A/2 at A > A, 
and ^ A in the opposite limit. It is important to observe 
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FIG. 1: (a) The f2g-orbital degeneracy of Co''+ (d^)-ion is 
lifted by trigonal distortion and spin-orbit interaction. A hole 
with pseudospin one-half resides on the /-level. Its wavefunc- 
tion contains both Eg and A\g states mixed up by spin-orbit 
coupling, (b) The hopping geometry on the triangular lattice 
of Co- ions. 0/3(7) on bonds should read as tap = t, t^-y = —t', 
and a,/?, 7 G {a,h,c} with a — dyz, b — dxz, c = d^y. The 
t-hopping stems from the charge-transfer process via oxygen 
ions, while t' stands for the direct d-d overlap. 
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that the pseudospin /g- states 

It)/ = zce| + l)U)-se|0)|T), 
II)/ = zc,|-l)|t)-se|0)U) 



(2) 



are coherent mixtures of different orbital and spin states, 
and this will have important consequences for the sym- 
metry of the intersite interactions, as we see below. 

We model the band motion by the following Hamilto- 
nian suggested by the edge-shared structure (171] : 



(3) 



where t = t^^/ IS. 



pd originates from the d-p-d process via 
the charge-transfer gap Ap^;, and > is the direct d- 
d hopping. On each bond, there are two orbitals (a, [3) 
active in d-p-d process, while the third one (7) is trans- 
fered via the direct d-d channel [Fig^b)]. The hopping 
geometry in real situation could, of course, be more com- 
plicated (e.g., include p-p hoppings). The t,t'-model is 
the simplest possibility chosen for illustrative purposes. 

The quasiparticle band structure is obtained within 
a slave boson approach, where the electron operator is 
represented as e\oiia => y/xctia+incoherent part. 

The last equation implies that we consider a fermi-liquid 
regime, where holons (described by the el operator) are 
condensed, while the cxia operators represent coherent 
fermionic quasiparticles. On the mean-field level (the in- 
coherent part being discarded), this leads to t tcoh — 
xt, and one obtains a renormalized hopping Hamiltonian 
jffcoh ^ xHf. While such a Gutzwiller-type approxima- 
tion fails to capture high-energy electronic processes far 
from the Fermi level, it is believed to provide a qualita- 
tively correct description of the low-energy states. 

Under the transformation the hopping Hamilto- 
nian obtains a rather complicated matrix structure in 
the full spin/orbital space, and the elements of this ma- 
trix sensitively depend on the angle 9. The Hamiltonian 
TJa + H\ -\- H^"^ has been diagonalized in momentum 



FIG. 2: (a) Quasiparticle dispersion curves in a slave-boson 
mean-field approximation, (b) Density of states (in units of 
1/t) near the Fermi level. The latter is indicated by thin gray 
lines, (c) Fermi surface in the Brillouin zone. V = (0,0), 
M = (0,27r/V3), K = (47r/3,0). In all the panels, solid 
(dashed) lines correspond to the doping level x = 0.2 {x — 
0.3). Parameters used: t' /t = 0.8, \/t = 0.7, A/t = 1.0. 



space numerically. The obtained quasiparticle dispersion 
curves are shown in Figl^^a). The first point to notice is 
that even for a substantial doping level, x — 0.3, there is a 
clear separation of bands which can be traced back to the 
on-site level structure discussed above. In particular, we 
find at small doping that the states near the Fermi level 
are derived dominantly from the /-pseudospin states; 
therefore, they are dispersive modes with mixed spin- 
orbital quantum numbers. The "/"-band has the bot- 
tom at the F-point [i^ , and shows flat portions near the 
if-points in the Brillouin zone of the triangular lattice, 
which lead to a singular density of quasiparticle states 
near the Fermi level [Fig[21[a,b)]. An interesting dop- 
ing evolution of the Fermi surface (FS) [Figl^fc)] is also 
related to the presence of those flat portions. The over- 
all energy window, covered by three quasiparticle bands, 
is reduced to ~ 3i by correlation effects. The high en- 
ergy, incoherent electronic states, extending up to a free- 
electron scale 9t, are beyond the present approxima- 
tion. The consideration of the incoherent part of the 
electronic motion, which is a difficult task in correlated 
models in general, is further complicated here because of 
the additional orbital degrees of freedom. 

As far as the low-energy physics at small doping is 
concerned, we arrived in fact at a single-band picture for 
states near the FS. However, this band has little to do 
with a conventional single-orbital band of Aig symmetry 
suggested by a free-electron band calculations ji^l and 
taken in Refs.Tol llllfl^llsl as a starting point to develop 
t-J model physics for cobaltates. The crucial difference 
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here is that the quasiparticle wave- functions are made 
of states which do not conserve neither spin nor orbital 
angular momentum separately; rather, they are build on 
eigenstates of the total angular momentum, as reflected 
in Eq.(|2I). This unusual situation results from the inter- 
play between strong spin-orbit coupling of Co*+ ion and 
the quasiparticle kinetic energy reduced by correlation 
effects. The main message is that spin-orbit coupling 
becomes increasingly effective near the Mott insulating 
limit and is thus essential for the formation of the quasi- 
particle bands, and hence for the FS topology. This im- 
plies that the FS shape at low doping may not necessarily 
be similar to that given by a free-electron picture, as one 
would expect in single-band systems like the cuprates. In 
a more general context, the outlined picture might be rel- 
evant also to other transition metal oxides, where narrow 
quasiparticle bands are derived from (quasi)degenerate 
t2g-states with strong spin-orbit coupling, such as com- 
pounds based on late-3d, 4d and 5d-ions. 

Now we turn to the superexchange interactions which, 
in analogy with high- Tc cuprates, could be one of the rel- 
evant interactions responsible for superconducting pair- 
ing. When the FS states are derived from pseudospin 
states with mixed spin and orbital quantum numbers, an 
important question arises about implications of such a 
mixing for the pairing symmetry. Projected on the pseu- 
dospin states, the superexchange interactions may in fact 
give nontrivial pairing channels, which were not present 
in the original, pure-spin Heisenberg model. 

A superexchange Hamiltonian in orbitally degenerate 
systems reads in general as Hj = J [{Si ■ Sj)Aij + Kij] 
|2J|. The energy scale is J = At^ /E, where the virtual 
charge excitation energy E is determined either by the 
on-site Coulomb Ud or the charge-transfer energy Ap^, 
depending on which one is lower. The operators Aij, Kij 
represent the orbital degrees of freedom. Without orbital 
degeneracy (e.g., in the cuprates) Aij = l,Kij = —1/4, 
and Hj supports uniquely singlet pairing. In the present 
situation, one has to: (i) derive a full structure of the 
orbital operators Aij and Kij (which usually depends on 
the (ij)-bond orientation in crystal via the hopping ge- 
ometry 13) I ^^'^ project the obtained Hj onto the 
active pseudospin /^-subspace given byEq.©. Details of 
this lengthy derivation ^] will be presented elsewhere; 
for our purpose, the result can conveniently be repre- 
sented in the following form: 

(4) 



Here, Jf = J/9 with J = 4t'^/Apd- In equation igj, 
the interactions are separated into pseudospin singlet and 
triplet channels, namely, Sij = {fi^fjl-Iilfj'[)/V^, while 



(5) 



with Uj^o = iiU^ffi + f,ifj^)/V2 , ^ ff^fj^ and 
tij^-i = filfjl- ^ phase if>s in Eq.© and below de- 
pends on the (ij)-bond direction S: (pg = (0, (j), —(f)) on 



5 = (12, 23, 13)-bonds (see Fig. lb), respectively. The rel- 
ative weights of the different components (M = 0, ±1 
projections of the total pseudospin of the pair) are con- 
trolled by the angle 9, defined above, via — {I + 
cos20)/2 and Uxy = sin26'. Finally, a constant 
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3 cos 261 - 1 



,(6) 



while for the triplet channel one has nt — 3/2. 

In general, the interaction |@J supports both singlet 
and triplet pairings on the fermi-surface. Certainly, there 
are also some other pairing forces in the cobaltates; yet 
it is interesting to explore the outcome of the superex- 
change interactions alone. We have therefore calculated 
the mean-field superconducting transition temperatures 
as function of the parameters involved in our effective t-J 
model for the "/"-pseudofermion band. 

In the singlet channel, it is known [lol Ull 
that pairing on the triangular lattice is optimized by a 
complex order parameter (sy ) = D e"^'' (the degenerate 
conjugate state is obtained by (fig — > —(f)g)- In momen- 
tum space, the gap function is determined by jd{k) — 
"ii{k) + ij2{k), where 71 — cos — cos ^ cos and 



This state breaks the time- 



72 = V3i 
reversal symmetry. 

As suggested by the very form of Eq.Q, the bond 
dependence of the triplet order parameter is best 
parametrized according to the projection AI of the 
Cooper-pair pseudospin. Namely, a positive interference 
among the different Af -channels is achieved by 



{tzj,±l) 



Ti-96 



(7) 



The pairing amplitudes dz , dxy are proportional to and 
axy^ respectively. The bond arrangement of the phases 
g-iM05 translates in momentum space into the form- 
factors ^z{k) and "f^y{k) = ±^x{k) + ijy{k), where 
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sm Kt — 2 sm — cos — 
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/— kx . ^/Sky 
V 3 cos — sm — - — , 

, , . kx ^/Sky 

sm kx + sm — cos 

2 2 



(8) 
(9) 
(10) 



The gap function has no nodes (apart from the F- and Ap- 
points), and the superconducting gap anisotropy is given 

by A{k) oc y |az7i:| "I" |'^a;i/7a:y| which depends on angle 
9. Remarkably, the pseudospin-triplet state is nondegen- 
erate and thus respects the time-reversal symmetry. This 
is because the orbital currents associated with M — +1 
and M = —1 components flow in opposite directions (ob- 
serve ± signs in above equations) , and the orbital angular 
momentum of the Cooper-pair is quenched. 
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FIG. 3: (a) Doping x and (b) trigonal field splitting A depen- 
dences of the mean-field superconducting transition tempera- 
ture Tc. Parameters used: t' /t = 0.8, \/t = 0.7, Jf/t = 0.05. 



physics between cobaltates and ruthenates |2g] and also 
the recently discovered superconductors KOS2O6 [29| . 
where Ru(4d)- and Os(5ci)-electrons with strong spin- 
orbit coupling reside on nearly degenerate t2g levels. 
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The mean-field Tc for the triplet state (assuming the 
holons are already condensed) is obtained from 



l = Jf 





1 ± 1^ 




\^xylxy\ 



k 



2CS 



tanh 



2r, 



(11) 



where is a "f '-quasiparticle energy. To obtain in the 
singlet channel, one simply has to replace the numerator 
in this equation by (^Hs\"fd{k)\^)- 

Numerical calculations show that for t' < t and 
^pd/Ud < 1 (relevant for cobalt oxides j2Q]) the triplet 
state is always favored. Shown in FigOl are some plots 
for Tc, calculated for an effective exchange parameter 
Jf = 0.05t (implying ~5 meV for t ~ 0.1 eV). Nonmono- 
tonic dependences of Tc on doping and the orbital split- 
ting (induced by c-axis compression) are due to sensitive 
variations of the quasiparticle band near the if-points of 
Brillouin zone (Fig|2J). We note that the obtained re- 
values (~ 5 K for t ^ 0.1 eV) and their dependences are 
qualitatively consistent with experimental data 0, . 

The superexchange driven triplet state is interesting 
on its own. This state is stabilized by spin-orbit coupling 
acting on the low-energy states of doped Mott insula- 
tor with t2g orbital degeneracy. Its basic features (no 
time-reversal symmetry breaking, only partial Khight- 
shift drop below Tc, high critical fields i?c2, etc) are 
reminicent to reported data in cobaltates. However, a 
specific comparison with (ongoing and still controver- 
sial) experiments would require that some other relevant 
physics (electron-phonon coupling, charge-ordering ten- 
dencies that are enhanced by orbital-polaron effects , 
phase separation, etc), may need to be included in a more 
realistic theory. We believe that the present work should 
be a proper starting point for such a theory. 

To conclude, the low-energy Fermi-surface states in the 
underdoped cobaltates are strongly influenced by spin- 
orbit interaction intrinsic to the t2g electrons of Co^+ 
ions. This leads to nontrivial symmetry of the superex- 
change interactions and to a novel pseudospin-triplet 
paired state, qualitatively different from that in the high- 
Tc cuprates. Rather, there might be some common 
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